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0^ . Carbon nanotubes provide a remarkably versatile system in which to explore the effects of 

O^ ' Coulomb interactions in one dimension. The most dramatic effects of strong electron-electron 

repulsion are ortbogonality catastrophes. These orthogonality catastrophes come in different 

(~| ■ varieties, and can be observed both in low-bias transport and tunneling measurements on 

^ ' nanotubes. This article contains a review of previous work and new material on the crossover 

between Coulomb blockade and Luttinger behavior. 

(N 

1 Introduction and Model 

j^H ■ In this section we describe the basic influence of interactions in carbon nanotubes, following 

1^ , Ref. 1. The band structure of metallic nanotubes has been discussed by several authorsB'S It 

■y I is well-captured by a simple tight-binding model of pz electrons on the honeycomb lattice. For 

the metallic tubes, evaluating the resulting tight-binding band structure for the discrete set of 
allowed quantized transverse momenta Qy leads to only two gapless one-dimensional metallic 
bands (for the simplest {N,N) armchair tubes, these have qy = 0)a'El These dominate the low- 
O ' energy physics, disperse with the same velocity, vp, and can be described by the simple Id free 

^,. Fermion model, 

^RiJdx^Ria - ^LiJdxi^Lia] > (1) 









Ho = y^l dxvF 
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CN I where i = 1,2 labels the two bands, and a =], [ the electron spin. We neglect curvature effects 

(which could open up small gaps in all but the armchair tubes), since these are subdominant to 
\^ . the Coulomb interactions. Higher sub-bands (for the armchair tube these have qy ^ 0) are of 

^P I course also present, and will be returned to in Sec. 3. 

We will make heavy use of the bosonized representation of Eq. ||, obtained by writing 
i^R/Liia ~ 6**^'^*"='=^'°), where the dual fields satisfy [(/>jq,(x), 6*^73(2/)] = —m5ij5a(iQ{x — y). Ex- 
J^ \ pressed in these variables (1) takes the form Hq = X^i « '^o(^icn i'ia) 



^ : n,{e,4>) = j dx'^Hdjf + {d^<t>f]. (2) 

^ ■ The slowly varying electronic density in a given channel is given by pia = ^ Ria'^ RioL~^''rLia^ Lia ~ 

dxOict/T^- The normal modes of TIq describe long wavelength particle-hole excitations which 
propagate with a dispersion uj = vpq- 
/\ ' Turning to the interactions, a tremendous simplification occurs when N is large: the only 

j^ . couplings which survive in this limit are forward scattering processes which involve small mo- 

mentum transfer. Roughly speaking, this can be understood as follows. "Interbranch" scat- 
tering processes (such as backscattering and umklapp) involve a momentum transfer of order 
Ikp ~ 1/a, where a is the carbon-carbon bond length. The matrix elements are therefore dom- 
inated by the short range part of the interaction, at distances r ^ a, where the interaction 
changes significantly from site to site. However, the electrons in the lowest sub-band are spread 
out around the circumference of the tube, and for large N the probability of two electrons to 
be near each other is of order 1/A^. For the Coulomb interaction, the resulting dimensionless 
interaction vertices are of order [e^ jhvp) x l/AoH. By contrast forward scattering processes, in 
which electrons stay in the same branch, involve small momentum exchange. They are domi- 
nated by the long range part of the Coulomb interaction, at distances larger than the radius, 
and there is no 1/A^ suppression. 



For A^ > 10 it is thus appropriate to consider a Luttinger model, in which only forward 
scattering vertices are included. A further simplification arises because the squared moduli of 
the electron wavefunctions in the two bands are identical and spin independent. All the forward- 
scattering vertices can thus be written as a single interaction, coupling to the total charge density 

Ptot = J2ia dxOia/'n-- 

We will suppose that the Coulomb interaction is externally screened on a scale Rg, which 
is long compared to the tube radius R, but short compared to the length of the tube. For 
simplicity, we model this by a metallic cylinder of radius Rg, placed around the nanotube. From 
elementary electrostatics, the energy to charge the nanotube with an electron density eptot is 

2 



Hi^t = e'HRs/R)J dxptf (3) 

Since H[^t only involves ptot it is convenient to introduce a spin and channel decomposition 
via, 0j,p/(T = (^jT ^ ^il)/V^ and 9^± = {di^ it 92f_i)/V^ with p, = p,a, and similar definitions 
for (p. As defined, the new fields 6a and (pa with a = {p/a,±), satisfy the same canonical 
commutators [4>a{x), db{y)] = —iTT6ab&{x — y). In the absence of interactions the Hamiltonian is 
simply Hq = J2a IxT'^o(^a,4'a)- which describes three "sectors" of neutral excitations and one 
charged excitation. Including the interactions only modifies the charge sector, which is described 
by the sum of two terms Tip = 'Ho(^p+) 4'p+) + -f^int(^p+) which may be written 



K, = / <l.^ 
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\dxep+f + g{dx^p+f] . (4) 



This describes the Id accoustic plasnion which propagates with the velocity 

Vp = \/vf{vf + (Se^/vrTi) ln{Rs/R)) and is characterized by the Luttinger parameter g = vp/vp. 

Taking a screening radius of Rg = 0.1/u, one estimates g ~ 0.2. 

2 Low-Bias Transport 

As evidenced in Eq. ^ the elementary excitations of the Luttinger liquid are collective modes, 
very different from the quasiparticles of a normal Fermi liquid. The most natural experimental 
measure of this non-Fermi-liquid physics are therefore those which probe the overlap between 
these excitations and bare electrons. The simplest way to achieve this is via a low-bias transport 
experiment in which electrons are removed from a metallic source lead and injected into the 
nanotube (and vice versa at the drain). Surprisingly, for this purpose it is beneficial to have 
poor contact between the leads and the tube. In this poorly-contacted limit, transport proceeds 
by sequential tunneling of individual electrons into and out of the nanotube, and is thus a 
measure of the single-particle overlap. By contrast, with ideal (adiabatic) contacts, the charge 
on the tube is no longer a good quantum number, and transport occurs via a collective fiow 
of charge. Despite the significant different between the single-particle spectrum of Luttinger 
and Fermi liquids, their collective mode spectra are quite similar. Indeed, as discussed by 
many authors, an ideally contacted pure Luttinger liquid has a two-terminal resistance which is 
completely unaffected by interactions. The ideal four-terminal resistance of a zero temperature 
Luttinger liquid is zero, a characteristic of ballistic transport. 

We consider the transport through a single tunnel junction, which could be between a normal 
lead and the nanotube, or between two nanotubes. The former situation has been investigated 
experimentally and theoretically in Ref. 6 (see the paper by Bockrath for details). In general, 
the current through a single tunnel junction can be determined perturbatively in the tunneling 
amplitude w. Consider the tunneling Hamiltonian 

HtunnXt) = E ^^''''4aCra + W* e-'^^'cUq^, (5) 



where (yj-Jc^i create/annihilate electrons to the left and right of the junction. If the right-hand- 
side is a nanotube, there is in general some tunneling amplitude into each of the Dirac channels: 

(r) 

Cra = J2i=i,2^p=R/L^Pi'^Pia- We have chosen to include the voltage difference between the 
two leads here as a time-dependent vector potential. The current across the junction follows 
from gauge invariance: / = dC/dax, where Ux is the component of the vector potential normal 
to the junction. The first non-zero contribution to the current in time-dependent perturbation 
theory in w is 

(/(t)) CC \W\^ f e^^(*-*')Aet.(t - t') ^ Aet.(^), (6) 



where ^ret. is the standard retarded correlation function of the operator A = cj^c^^. Since 
the overall magnitude of the current depends on an unknown transmission probability, we will 
simplify various formulae by omitting constant prefactors where possible. The retarded corre- 
lator is obtained in the usual way by analytic continuation from the imaginary-time correlator 
A{t) = {TrA{T)A^O)). In particular, 

dte-'^U{^-it), (7) 

where, as usual (3 = l/ksT. 

2.1 Thermodynamic Lim,it 

We first consider the case in which the systems on either side of the contact are semi-infinite. 
This is a valid assumption for a finite Luttinger liquid provided the temperature or bias voltage is 
large compared to its capacitive charging energy. For both the case of Fermi-liquid to nanotube 
and nanotube-nanotube junctions, the imaginary-time correlator above takes the form 



_ N^ 2+a 

TTT 



with a > 0, which can be obtained straightforwardly from the bosonization formulae in Sec. 1. 
This form obtains in many distinct cases - e.g. bulk contact between a Fermi liquid and a 
nanotube, bulk contact between two nanotubes, and end-to-end contact between two nanotubes. 
The different cases are distinguished by the value of the exponent a. To get a feeling for the 
physical meaning of a, imagine a Fermi's golden rule estimate of the current across the junction. 
If the initial state is the ground state, the final state consists of e.g. a hole in the left lead and 
an electron in the right. Fermi's golden rule gives a density of states factor which is the product 
of the electron density of states on the right and the hole density of states on the left. Indeed, a 
direct calculation of the tunneling density of states p*')"' for a Luttinger liquid gives the formula 

a = ai + Or, where p*T'(e) ~ |e|"''/'- 

The density of states exponent depends crucially on the Luttinger parameter g, and also on 
the point of contact. In particular. 



(g ^ — l)/4 near cap 

(ff + 5""^-2)/8 in bulk 



«=U! "t^^o^.o ::tr" (9) 



Note that, since g < 1, the orthogonality exponent a is significantly larger near the end of a 
nanotube. This enhanced interaction effect arises essentially from the decreased ability of an 
added charge to spread away from the cap area relative to the bulk. A contact may be considered 
an end contact if the distance from the end of the nanotube is smaller than the lesser of the 
characteristic length scales Lt = Hvp/kBT, Ly = hvp/eV. It is important to realize that this 



is not a small effect: for g <^ 1, the orthogonality exponent near the cap is nearly twice that in 
the bulk. 

Using Eq. g and Eq. 0, one obtains the I-V curve 



/ = loT'-^" sinh(/3y/2) 
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(10) 



where Iq oc \'w\'^ is an unknown prefactor and T{z) is the Gamma function. This form is 
well-known from the literature on tunneling and Caldeira-Leggett models. A slightly more com- 
plicated formula can be obtained for the differential conductance G = dl/dV by differentiation; 
differentiating only the smh{fiV/2) prefactor gives a common approximate result (as the remain- 



der is small at low biases). Eq. 10 has an important scaling property. The quantities J/T^+°' 
and G/T" are functions only of the ratio V/T. This implies that IV curves taken at different 
temperatures should collapse if, e.g. G/T"' is plotted as a function of V/T. See the paper by 
Marc Bockrath in this volume for an experimental observation of this behavior. 



2.2 Charging Effects 

Coulomb blockade is a ubiquitous phenomena in mesoscopics physics. For a few micron- long 
nanotube, the capacitive charging energy can be of the order of tens of Kelvin, and affects 
the transport dramatically at low temperatures. The crossover between Luttinger liquid and 
Coulomb blockade behavior is a complex and relatively unstudied problem. We will discuss 
a simple model of charging effects for a single tunnel junction between two nanotubes. The 
charging behavior is, unfortunately, much less universal that the Luttinger liquid results of the 
previous subsection, so we must make a number of assumptions to progress. 

A panoply of energy scales exist in a finite Luttinger liquid. The largest is generally the 
charging energy. In the Luttinger model, this is the energy of a zero mode in which each of the 
phase fields Oa winds by ^ti 12. For g' <C 1, this is dominated by the charge mode ^p, and from 
Eqs. ^^ Ec ~ nhvp/Sg'^L ~ {e^ /L) In Rg/R. For a one-dimensional system, the level spacing 
also scales inversely with the length. In particular, the Luttinger liquid exhibits two such level 
spacings due to spin-charge separation: a plasmon energy Ep = Trhvp/L and a single-particle 
energy eq = -rrhvp/L. For g <^ 1, one has Eq ':$> Ep ^ Eq. We will therefore consider only the 
largest of these energies, Eq- This corresponds physically to including charging effects but not 
individual level quantization. 

Within this charging-only model, we will consider the specific case of a single tunnel-junction 
between two finite-length nanotubes. The treatment here is valid provided that the resistance 
of the nanotube-nanotube junction is much larger than the resistance to the leads. In this case, 
at least for not too low temperatures, transport is dominated by the internal tunnel junction. 
We require i?junction » -Rioads » h/e'^. 

To proceed, we must specify the capacitances in the system. In general, the charging energy 
is given by an energy function Em ,n2 1 where ni and 77-2 are the charges (in units of e) of the left 
and right sides of the system. Given such a function, Eq. |^ can still be applied, but A{t) must 
be recalculated. Because the charges are zero-mode quantities, the statistical average simply 
factors into a product of a zero-mode contribution and the thermodynamic contribution in Eq. |8|. 
Inserting into Eqs. &|^, one obtains 



/ = /o^— ^sinh(/5F/2) J2 e-f^^"i'"2 

ni,n2 



1 + - + ^^ 



AK 



ni,n2; 
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(11) 



where the "partition function" is Z 

^il,n2 = 2 (-E'ni+l,n2-l + E^^ni), ^E, 



ni,n2 



J2nin2^ /3^ni,n2^ "YYg have also defined the energies 

In general, a quadratic energy 



E., 



■ni+l,n2-l " 



-E, 



ni,n2- 




Figure 1: Ground state charge configurations in the n-rn plane, assuming W <U . Inside each hexagon the charge 
state is indicated by the ordered pair (711,712). For W > U, the pattern is rotated 90 degrees and there are no 
horizontal segments. Some possible trajectories as the gate voltage is varied are indicated schematically by the 

dashed red lines. 



function £'ni,7i2 contains five non-trivial parameters. It is helpful to consider a simplified form, 



E, 



7^\2 



711,712 



U{ni + 77-2 — n) + W{ni — n2 — m) 



(12) 



Here W can be viewed as a capacitive charging energy for the junction itself, while U gives an 
"bulk" contribution. A gate potential is included via n,m. Changing a single gate potential 
corresponds to motion along some straight line in the n^m plane. Eqs. 11-12 encompass a 
tremendously rich behavior. Consider the Coulomb blockade regime at low temperature. At 
zero temperature, the ground state charge configurations of the system (for w = Q) are shown in 
Fig. 1. As m,n are varied, the system makes two types of transitions between different charge 
states. Across the diagonal lines, one electron is added or removed from either the left or right 
nanotube, changing the net charge (and spin) on both tubes. Across the horizontal lines, one 
electron is transferred from the right to the left tube or vice versa, leaving the total change 
unchanged. In the limit we are considering, transport is dominated by tunneling across the 
inter-tube junction. Peaks in the zero-bias conductance thus occur only across the horizontal 
lines. 

As an external gate potential is varied, the system traces out a straight line in the m-n 
plane. The particular slope and intercept of this line is determined by the precise geometry of 
the tubes and the gate. The system will thus generically undergo a sequence of ground-state 
transitions, some of which (crossing the diagonal lines) do not give rise to conductance peaks, 
but which change the charge and spin. Such "internal" transitions have indeed been observed 
in nanotubesQ Moreover, depending upon the particular path taken, the locations of Coulomb 
blockade peaks as a function of gate voltage (crossings of horizontal lines) can be quite complex, 
e.g. quasiperiodic or containing a fairly regular series of peaks followed by a region of zero 
conductance. 

Let us return to the problem of the crossover between Coulomb blockade and Luttinger 
behavior, contained in Eqs. |ri|-12. A general analysis is beyond the scope of these proceedings. 
We therefore content ourselves with a few simple situations. The most "universal" limit is the 
high-temperature case, basically when T > W (note that W and not U sets the size of the energy 



that an electron must obtain to cross the tunnel junction). In this case, the sum in Eq. 11 can 



be approximated by an integral, and 

At very high temperatures, the exponential terms in W are negligible and the result reduces to 



Eq. 10. An appealing feature of Eq. 13 is its independence of m,n and U. 



Unfortunately, the low-temperature limit is much less universal. As an illustration of the 



application of Eqs. 11 12, consider the behavior of the zero-bias conductance G(T) as a function 
of temperature when the gate voltage is tuned to a zero-bias conductance peak. In this case, 
one finds power-law behavior at both low and high temperatures: G{T) ~ Go{W/T)T'^. The 
prefactor, however, has the amusing property that Go{W/T » 1) = Go(W/T <C l)/2. The 
factor of two comes physically from the fluctuations of the system between two generate ground 
states at low temperatures, and reflects the zero-temperature entropy of the two degenerate 
ground states of the system! Amusingly, a more detailed examination shows that G{T)/T°' 
actually crosses non-monotonically between these limitsp 

3 High-Bias Tunneling Spectra 

We have seen in the previous section how Coulomb interactions give rise to numerous interesting 
non-linearities in low-bias tunneling transport. Tunneling spectra also provide a powerful probe 
of high-energy excitations. In nanotubes, tunneling experiments have verified the presence 
van Hove singularities in the density of states (at eV energies) due to the higher sub-bandsBI 
In non-interacting band theory, these singularities are of inverse square-root form, giving a 
contribution /9o(e) ~ \/rn{e — A)~^'^Q{e — A) for energies just above the subband edge at e = A 
{m is the subband effective mass). 

How do interactions affect these van Hove singularities? A simplified, though unphysical 
model in which the mass of the higher subband is taken to be infinite provides considerable 
insight. In this limit the higher energy subbands can be replaced by discrete, localized levels. The 
"x-ray edge" problem of a localized level interacting with a conduction sea was solved by Nozieres 
and de Dominiciill, and is one of the first demonstrations of an orthogonality catastrophe. 
Physically, the core hole is "dressed" through interactions with conduction electrons, which see 
the hole as a scattering center. This leads to a broadening and reduction of the tunneling density 
of states from a sharp delta-function to a power law singularity. 

RemarkabljjJ^Jiese x-ray edge effects persist even for this finite mass case, as we now proceed 
to demonstratep3'EJ We argue that a necessary and sufficient condition for the presence of such 
finite-energy singularities is a conserved quantum number distinguishing the states of the higher 
subband from the conduction states. In the case of the carbon nanotube, this is an angular 
momentum quanta. If such a distinguishing quantum number is absent (as might occur in a 
nanotube due to breaking of the rotational symmetry by interactions with a substrate), we 
expect the van Hove peak to be rounded and rendered completely nonsingular. 

We describe here a simple forward-scattering model of the interaction of the conduction 
electrons with the higher subband - other interaction channels are discussed in Ref. 13. Near the 
putative van Hove singularity, the unoccupied Id subband can be described by a non-relativistic 
electron operator d,d': 

\^. (14) 



<= jdxdi 



-l^dl + l^ 
2m 



Here A is the gap to the first subband and m is an effective mass. The electron field satisfies 
{da{x),dp{x')} = 5a(i6{x — x'). In the case of a carbon nanotube, there are actually multiple 
degenerate subbands at energy A. This degeneracy is unimportant within the forward-scattering 
model, as the tunneling DOS involves only states with a single excited electron. 



The interaction Hamiltonian is ifint = J^'Hmti with 



Hint = e'ln(i?./i?) 



(d)df + -dJp+Sd 



To understand the effects of Wint. , consider the canonical transformation 
ep+{x) = ep+{x)--f f dx'd\x')d{x'), 



dix) 



Ji4>p+(x)/t^ 



d{x), 



(15) 

(16) 
(17) 



where 7 = 7r(l — g'^)/2. Eqs. 16-|T7 embody the physical process in which the conduction sea 
adiabatically adjusts to the heavy particle. In particular, Eq. 16 represents the depletion of the 



conduction electron density near the heavy particle due to Coulomb repulsion. Eq. |l^ represents 
phase shifts of these conduction electrons when the heavy particle is introduced. Formally, the 



exponential of the dual (</>) field in Eq. 17 is a Jordan- Wigner "string" operator which has been 
attached to the heavy particle. 



Remarkably, although the canonical transformation, Eqs. 16-17, does not remove the inter- 
action completely, it does transfoim the Hamiltonian into one with only irrelevant couplings in 
the renormalization group senseEj This indicates that at long times and distances, the trans- 
formed fermion and boson correlation functions asymptotically factorize. From this result, one 
straightforwardly obtains a modified van Hove singularity, gives a modified van Hove singularity: 



Pie) 



Po 



A 



'P 



e(e-A), 



(18) 



where Q{x) is the heavyside step function, and the orthogonality exponent j3 = {1 — g'^)^/{8g) ~ 
0.3 for typical metallic nanotubes. 
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